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Abstract. We investigate the coupling of an inhomogeneous electron system to 
phonons. The properties of an electronic system composed of a mixture of microscopic 
ordered and disordered islands are changed fundamentally by a phonon mode. In 
high-Tc cuprates, such a phase separation scenario is supported by recent STM and 
nuclear quadrupole resonance studies. We show that even a weak or moderate electron- 
phonon coupling can be sufficient to produce dramatic changes in the electronic state 
of the inhomogeneous electron gas. The spectral properties calculated in our approach 
provide a natural explanation of the low-energy nodal ARPES features and exhibit a 
novel non-Fermi-liquid state stabilized through electron-phonon coupling. 
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1. Introduction 

An eminent feature observed in the angle-resolved photoemission spectroscopy (ARPES) 
studies of high-Tc cuprates is the kink-like change of electron velocity, interpreted in 
terms of the coupling to oxygen phonon modes [H |2] . These kinks are associated with 
a distinct break-up of the spectral weight into a high-intensity part which develops 
near the Fermi surface, and a broad structure of less intensity at higher energies below 
the Fermi level. Despite the structural differences between various types of cuprates, 
the peculiarities in the dispersion at 50-80 meV were detected in Bi2Sr2CaCu208+5 
(Bi2212), La2-,;Sr^.Cu04 (LSCO), YBa2Cu306+x (YBCO) and other related systems 
[H [21 m [5l m [6] , and are believed to shed light on the microscopic mechanism of high-Tc 
superconductivity. Although significant efforts have been directed towards a theoretical 
investigation of the influence of many-body interactions on the one-electron properties 
of Hubbard and t-J models, typically used for cuprates [7J, the origin of the nodal 
ARPES features and of their unusual doping and temperature behavior still remains an 
open question. 

Many-body effects are known to produce changes in the electronic dispersion, as 
deduced from electronic photoemission spectra. One mechanism, which was shown to 
cause low-energy kinks, is controlled by a strong electron-phonon coupling in the free 
electron gas [HI [9]. This coupling leads to a renormalization of the electron effective 
mass in the energy range close to the Fermi energy Ep, below a characteristic phonon 
frequency tUph: E — Ep < cjph. On the other hand, when the electrons in a metal are 
strongly correlated, a purely electronic mechanism can also lead to dispersion kinks 
which have been related to a crossover between Fermi- and non-Fermi-liquid behavior 

m- 

In the analysis of ARPES intensities, one should always consider the fact that 
photoemission experiments analyze electronic surface and subsurface states which brings 
the near-surface correlations into the focus. Due to the strong interaction of these 
electrons with subsurface phonon modes and virtual charge transfer excitations, the 
effect of such collective modes on the electronic subsystem is crucial [HI [12]. A direct 
consequence of these interactions is a significant reduction of the local Hubbard repulsion 
to values well below the electronic bandwidth 8t, where t is the electron hopping. This 
finding allows to suggest that mechanisms different from the purely electronic could 
play a role in the appearance of the dispersion kinks in cuprates. As the coupling with 
magnons cannot satisfactorily explain the doping and temperature behavior of such 
kinks, we will focus on the analysis of the interaction with phonons. 

In the high-Tc cuprates another important factor, which has to be addressed 
in studies of their electronic properties, is charge inhomogeneity. In the pseudogap 
state and in regimes with suppressed superconductivity, scanning tunnehng microscopy 
(STM) experiments indicate a local electron order p^[Til[T5l[T6l[2n[22l[23l[2l]. 
This, together with nuclear quadrupole resonance (NQR) and resonant soft X-ray 
scattering studies [25l |26l [271 [2H], provides strong support for a state of electronic 



Electron-phonon coupling in a two-dimensional inhomogeneous electron gas 3 

phase separation as one of the widely discussed scenarios for under- and optimal-doped 
cuprates [131 [HI [13 [251 l26l [27] . Such a state exhibits a mixture of microscopic charge 
ordered (characterized as "more insulating" with suppressed local density of states) and 
charge disordered uniform (metallic) domains with a dominance of the metallic phase 
at higher doping levels [Ej- The ARPES intensities obtained for the electronic phase 
separated state inevitably contain the contributions from both types of domains. This 
fact can be easily understood since the electrons collected on the ARPES analyzer can 
be emitted from the disordered as well as from the ordered surface islands. Due to 
the quantum nature of the collected electrons, they cannot be described in terms of 
pure "ordered" or "disordered" electron wave functions, but rather as a superposition 
of both states. Therefore, in order to understand ARPES intensities within such an 
inhomogeneous scenario, we need to analyze the consequences of the electron-phonon 
coupling not only for the disordered, but also for the ordered sections of the surface. 

In the hole-doped cuprates, the recent ab-initio studies have demonstrated a 
weakness of the electron-phonon coupling which should result in a negligibly small 
contribution to the formation of the dispersion kinks [T71 [18] . Despite the different 
approaches, these works have been focused on the electronic homogeneous state, without 
any serious attention to a possible electronic charge order. As a consequence of such 
a predomimant consideration of the uniform electronic state, the general view about 
a relative unimportance of the electron-phonon coupling for the electronic properties 
of the copper oxide planes became widely accepted in the literature. In an alterating 
approach, the strongly inhomogeneous interactions of the electrons to several optical 
phonon modes (buckling and breathing) has been proposed [19] which was not really 
successful in the explanation of real physical mechanisms of spectral anomalies in the 
cuprates. 

In the present work, we analyze the electron-phonon coupling in the charge-ordered 
state. We show that in contrast to the uniform electron gas, in the ordered system even 
a weak or moderate coupling to an optical phonon mode produces dramatic changes in 
the electronic properties and leads to a formation of new electronic state which cannot 
be described by a standard Fermi-liquid theory. Due to significant advances achieved 
in the recent STM studies of the cuprates [HI [20] , the existence of local electron order 
in these systems becomes rather a well established fact. To make a step towards an 
interpretation of the spectral features in the context of the inhomogeneous state observed 
in the STM studies, the inhomogeneous surface in our work is described as a mixture 
of the electronically uniform and ordered states. We obtain that the superposition of 
these states produces several distinct features and a characteristic intensity break-up in 
the electronic spectral maps. In our work, we connect this break-up with the spectral 
anomalies observed in the cuprates. Consequently, the proposed approach is expected 
to shed new light on the mechanisms of kinks and on the origin of non Fermi-liquid 
behavior detected in these systems. 
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2. Electronic spectral properties in the charge-ordered state 

To gain deeper insight into the ordered surface state, we provide a comparison of a 
disordered two-dimensional electron gas, characterized by a free tight-binding dispersion 
eo{k), with an ordered electron system, both coupled to an optical phonon mode with 
a frequency cUph = ujq. In our studies, we take into account the Coulomb interaction V 
between nearest neighbors 

Hei = T.k,a£o{k)nk^ + Vll(ij)ninj (1) 

where are the electron number operators and Ui = T^fjUifj. The electron-phonon 
interaction is considered in terms of a Holstein approach 

ifei-ph = -gT.ini{bi + bl) + UoT^iblk, (2) 

where the phonon operators refer to a vibration mode of frequency uq- The parameter 
g = ^UoEp {Ep is the polaron binding energy) refers the coupling of the holes in 
the copper oxide planes to the motion of apical oxygens in the top surface planes of 
the samples. In our calculations, the low-energy phonon frequency ujq = 0.05 — O.lt 
reflects the softening of the surface optical phonon modes suggested in Ref. |12j . 
For the free electronic dispersion eo{k) = —2trj^ — ^k (where r^^ = cosk^ + cos ky, 
= + 4:t' cos k^ cos ky and fi is the chemical potential) we choose t = 0.18 eV, 
t' = —OAt, which is in the range of typical values found from fitting ARPES data for 
Bi2212 and T12201 |29l|30]. 




Figure 1. Schematic view of checkerboard ording in a two-dimensional electron gas. 

In our analysis we consider a checkerboard electronic ordering which introduces a 
doubling of the square unit cell shown in Fig. [H This ordering is parameterized in terms 
of sublattice electronic occupancies n± with order parameter 6 = — and average 
electron concentration (n+ + n_)/2 = 1 — x. The values of 6 are obtained from the 
minimization of the mean-field type free energy. In such a charge ordered state, the 
electronic band structure is characterized by two subbands ej{k) = — -|- AVn ± A^. 
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The correction As = y {2V6)^ + {2tr]^ — Tj^^Y originates from intersub lattice electron 
correlations, introduced through the Hatree-Fock self-energy T,^^ = — J dq J^uj^ ^ " 
Vk-q ' G^^{q,ujn) in the Matsubara electron propagators G^^{q,u!n)- The new feature 
arising in the ordered band structure is the flatness of the energy subbands Bj{k) which 
is generated by their splitting through the charge-order gap Aq = 4:V6, Fig. [2]^a). Here, 
in contrast to the disordered free dispersion eo{k), the emergence of the gap Aq results 
in the formation of new local extrema of the ordered dispersions Sj{k). Fig. [3] shows a 
comparison of the detailed {k^., ky) map of the low-energy ordered subband e2{k) with 
the corresponding map of the free electron dispersion eo{k). One can immediately see 
that the electron order produces dramatic changes in the topology of the electronic 
structure which include (i) the formation of a new maximum in the central nodal 
point N = (71/2, 71/2) of the Brillouin zone; (ii) the appearance of new saddle points 
S = {k^,ky). The saddle point S is located close to F (the second symmetric saddle 
point is close to Y) in the nodal direction so that k^ = ky = ks where the coordinate 
k^ is given by 

^{At + 2VU^Y-{m'V5Y 

= — 8\t'm+2vi,.) — ' 

where the quantity parametrizes the Hartree-Fock self-energy T,'^^{k) = 

— (V/2)?7^J+_. For the chosen values of the model parameters, the maximal values 
of /+_ are of the order 0.5. The expression ([3]) clearly shows that the saddle point S 
is controlled by the charge order gap and disappears for weaker Coulomb repulsion V, 
i. e., it is absent in the disordered uniform state with 6 —>■ 0. We also note that the 
nodal point N indicated in the top map of the band eo{k) in Fig. [3] corresponds in fact 
to the intersection point of the two branches of £o(^) (see Fig. [2](a)). They are generate 
with the convolution to the first Brillouin zone and do not show an extremum at N. In 
contrast to eo{k), the new local extrema N and S of e2{k) are well-defined and should 
be considered as a generic feature of the band structure in the electron- ordered state. 

When the fiat low-energy subband 52 (fc) is coupled to the phonon mode, the local 
maximum at N and the saddle point at S lead to van- Hove singularities in the sublattice 
contributions to the electronic self-energies = —g'^fdqY,uj^,G^^{k — q,Un — 
Un')D{q^n')- Here D{q,uJn') = 2uJo/ {{iuJn'Y — uJq) is the unperturbed phonon Green 
functioiiji and the one-electron sublattice propagators G^^ {a = {+,—}) are calculated 
in the ordered electron state in the self-consistent Hartree-Fock approximation. 
Furthermore, the phonon scattering term can be decomposed as = Eq''=i=ASo'^ where 
the part AEq'^ ~ VS disappears in the charge disordered state. With this form of S^f", 
the renormalized electronic propagators can be conveniently presented as a combination 
of the band Green functions gj = l/{iun + ik — ^j{k,u!n))- Gaa = Co =F ^Gq, where 
Co = ((71 + 5'2)/2 and AGo = {2V6 + /S.T]^){gi — g2)/{T,i — ll2)- In the band propagators 

X As the main results for the electronic properties, including the topological anomalies of the dispersion 
and the NFL-state, appear already at small electron-phonon coupling Ep/t ^ 0.2, undressed phonon 
Green functions have been used in the calculations of the electronic self-energy. 
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gj, the effective self-energy parts S^- = + 4Vn ± ^(AS|;^ + 2V6y + {2t7]^ - S^^)^ 
contain the information about both, charge density fluctuations and electron-phonon 
scattering processes. We note that S2 and g2, which determine the low-energy quasi- 
particle excitations, are of prime importance. 

In Fig. El^b), we show the frequency-dependent real part of Sq'^(A;,C(j) at different 
doping levels x. The kinks appearing in T^q^ at u = Ui {i = 1, . . . ,3) correspond to 
the van-Hove singularities in S^f . In the charge disordered state, the typical van-Hove 
singularities in Ti^^ are caused by the local extrema in the electron dispersion eo{k) which 
appear in the high-symmetry points F, Y and M of the Brilloiun zone (Fig. [2]^a)). In the 
charge ordered state, the flattening of e2{k) produces the additional nodal maximum 
of e2{k = N) and the saddle point at = S which lead to the appearance of new 
van-Hove singularities in at small u. These low-energy singularities correspond 
to the peaks in ImSg'' (see Fig. Hj) and are related to the kinks in ReSf (Fig. [^b)). 
In the ordered state, the energies of the low-frequency singularities are determined 
by the equations u — £2(^4) ± = where ki = {r,N, S}. Specifically, while the 
maximum at = N produces a jump of Eg'' described by ImEg'^ ~ 6(fi'~^ — uj) at small 

= e2{k = N) + ujq > 0, the saddle point S leads to a distinct logarithmic singular 
behavior of T,^^: 

ImSg'^ ~ g^l + 60 - /(^ - ^0)) - Lu)/t\ 

+ ibo + fiu + ujo))\ogml-uj)/t\], (4) 

where bo = b{uJo) and /(cu) are the Bose and Fermi distributions functions, and the 
energies Q2/3 = ^2{ks) T ^0 are located close to the Fermi level: < fif < UJ2, 
UJ2 < 0. 

As the chemical potential directly affects £2 in these equations, the kinks of S2(fc, uj) 
at a; = Ui which result from the new van-Hove singularities are strongly doping- 
dependent. It is noteworthy that these low-frequency kinks at u = Ui {i = 1, . . . ,3) 
shown in Fig. EJ^b) will inevitably change the one-electron spectral properties near the 
Fermi level. To see this, we present in Fig. [2](c) the low-energy dispersion E2{k) of the 
underdoped system (x = 0.11), now calculated from the equation u + ^k — S2(fc, ou) = 
for the poles of the renormalized electron propagators for different values of electron- 
phonon coupling. A central property resulting from the new van-Hove singularities is a 
topological reconstruction of E2{k) close to the N-point. In Fig. [2]^c), (case Ep = 0.8t), 
this reconstruction corresponds to a transformation of the maximum E2{k = N) into a 
new singular hat-shaped band region which appears above a critical value E* ^ 0.2t. 

The novel nodal dispersion exhibits a nonanalytic character in the symmetric 
regions ni and n2 near the Fermi- level. In these regions, E2{k) is a multi- valued 
function of k. The two extra subbranches of E2 emerge from the additional poles 
of the one-electron Green function and are caused by the van-Hove kinks in Eo''(fc,a;). 
In the singular points 5*1 and S2 of the regions rii and n2 which connect monotonically 
decreasing and increasing branches of E2{k), the nodal electron velocity v = \'VkE2{k)\ 
approaches infinite values. The range k close to kp with singular v is separated from 
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Figure 2. Electronic band structure and self energy for a charge ordered system 
coupled to phonons where V/t = 1.3. (a) Band structure with vanishing phonon 
coupling in the ordered and charge disordered uniform states, where the dotted lines 
refer to ±uJo ^ 0.5t. Here kT/t = 0.03 and x = 0.11. (b) Real part of EP''(fc,w) 
calculated in the ordered state with ujq = 0.05t for different doping levels. Here k is 
located in the nodal region close to the point N = (7r/2,7r/2) of the Brillouin zone; 
fc| = 0.9, kT/t = 0.03 and Ep = 1.2t. (c) Evolution of the rcnormalizcd low-energy 
band structure in the ordered state with increasing Ep for x = 0.11 and loq — 0.05i. 
The transition at Ep = E* ^ 0.2t is characterized by a transformation of the local 
maximum E2{k — N) into a hat-shaped structure and by the buildup of a singular 
i?2(fc) in the regions ni and n2 with singular electron velocity v. The regions k ^ kp 
are separated by a gap A^r from fe N with hole- like excitations. The Fermi level is 
indicated by dots. 



the region A; N (characterized by hole-hke excitations) through a gap Ajv/t ~ 0.2. 
These topological anomalies are especially significant in the under- and optimal-doped 
range and for polaron energies above E*. The nodal properties of the new electronic 
state, stabilized at Ep = E*, cannot be classified in terms of Fermi liquid theory. 

In our analysis, the nodal non-Fermi liquid (NFL) behavior is also evidenced 
by an anomalously high scattering rate Ffc^ = Zfc^|ImSo''| (Zfc^ is the calculated 
quasiparticle residue). In Fermi- liquid theory of the free electron gas, one always finds 
llmEg'*! -C E2{k) sufficiently close to the Fermi level which signifies the existence of 
coherent quasiparticles with a long life time r = l/Fj,. At low T and very close 
to the Fermi level, the interactions of the free electrons to an Einstein phonon lead 
to an exponential decrease of ImEg'' ~ exp(— /^c^o) as T — 0. This means that 
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Figure 3. Electronic contour plots of eo{kx,ky) and e2{kx,ky) in the disordered 
{V/t = 0.2) and ordered (V/t — 1.3) state where the local extremal points are indicated 
by r, M, Y, N, and S. Here kT/t = 0.03, x = 0.11, Ep/t = 1.2. 

for small excitation energies e close to the Fermi level, there always exists a low- 
temperature range for which llmEg'^l <C e, which is consistent with the concept of 
Landau quasiparticles [31] . In contrast to this, in the ordered system the nodal 
electron-phonon scattering near the Fermi level produces a substantial ImEg'^ ~ 
5'^{log |fif /t| +log |^^f/t| +'r]N{t, t', V)} exp{—(3uJo) where the function rj^it, t\ V) results 
from a quadratic expansion of e2{k) in the vicinity of the nodal van-Hove singularities. 
Due to the large dominant contributions of these van-Hove singularities, in the 
temperature range kT/t < 0.03 which is of relevance for the cuprates, we obtain 
|ImSo''| > e where the excitation energies e < kT are located close to the Fermi level. 
The comparison of the corresponding ImEg'' in the free and in the ordered electron 
gas is presented in Fig. HI As a consequence, the strong increase of ImSg'^ and of 
the electron scattering rate Tj^^ leads to a violation of the condition for well-defined 
quasiparticles in the electron-ordered system. We emphasize that the ultimate reason 
for such a high electron-phonon scattering rate is the anomalously high density of the 
nodal electronic states related to the flat electron dispersion. The incoherent NFL nodal 
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Figure 4. (a) Imaginary part of Eq (fe, oj) and (b) scattering rate = llm^o I 
{Zkp is the quasiparticle residue) in the charge ordered state (case V/t = 1.3, 
|fc| = 1.11) and for the free electron gas {V/t = 0.2, |fe| = 1.3). Here k is located 
in the nodal direction near the N-point of the Brillouin zone, kT/t = 0.03, x = 0.11, 
and Ep/t = 0.2. 



state which develops for Ep > E* is in striking contrast to the Fermi-hquid state with its 
distinctive hole-type quasiparticle excitations for Ep < E* (dispersion similar to e2{k) in 
Fig. [2](a)). The existence of the singular multipole structure of the one-electron Green 
functions and of the anomalously high scattering rates has been confirmed by more 
elaborate calculations which include higher-order vertex corrections into the electron- 
phonon contribution to the electronic self-energy. As follows from these calculations, 
the inclusion of vertex corrections leads to slight shifts of the low-energy van-Hove 
singularities in ItoEq^, without qualitative changes in the nature of the NFL state. In 
fact, evidence for a NFL behavior along the nodal region, which is even unaffected 
by the onset of superconductivity, has been found for optimally doped BSCCO [1], in 
agreement with our findings. 

3. Electronic phase separation state 

As a state of electronic phase separation would imply a superposition of the ordered and 
disordered electronic states on the analyzer, the resulting ARPES intensity contains a 
combination of the disordered {A'^{k,u)) and the ordered {A°{k,u))) spectral functions: 
I{k,uj) = {cordA"{k,u) + CdisA°'(fc, a;))/(ti;). Here the coefficients Cord and Cdis refer 
to the ordered and disordered surface fractions. We assume that the ratio Cdis/cord 
depends linearly on doping so that Cdis/cord ~ x. The consequential dominance of the 
charge disordered contribution for larger x is consistent with the expansion of metallic 
disordered regions at higher doping levels as observed in Ref. [T6] . 

Fig. [5] shows the spectral intensity map in the nodal direction calculated for x = 
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Figure 5. Spectral intensity map I{k,uj) for an inhomogeneous electron system with 
X — 0.06 along the nodal direction. The bottom panel displays the detailed structure 
of the break-up feature at « —50 meV Here the spectral intensity is renormalized 
by the weight coefficient of the ordered state Cord; V/t — 1.3, kT/t = 0.03, loq — 0.05t 
and Ep — 1.2t. The position of the Fermi level is indicated by dashed line. 
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dis ord 



x=0.11 
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dis ord 
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0.2 0.4 0.6 0.8 1.0 1.2 

Momentum |k|/7i (nodal direction) 



Figure 6. Evolution of the spectral intensity map I{k,uj) at different doping levels. 
The position of the Fermi level is indicated by dashed line. 
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0.06. One can clearly distinguish two structures on this map. The needle shaped high 
intensity structure (1) is related to the free electronic band Eq. Moreover, a symmetric 
broad feature (2) appears, originating from the ordered state. The structure (2) is 
located on the plot in the range between ~ —60 meV and ^ —10 meV. The 
high spectral intensity of the structure (2) appears due to the singular character of 
E2{k) in the vicinity of the point Si of the nodal ni-region in Fig. [2t^c). Near the Fermi 
level, the NFL -region ni forms due to an additional van-Hove singularity at — lu^ 
(see Fig. EJ^b)) which results in multiple poles of the electronic propagator g2 and in 
high spectral intensity in the region between and uj°. In the intensity I{k,uj), the 
ordering-induced feature (2) is superimposed with the free band (1) which produces a 
break-up of the intensity at u; ~ —50 meV, shown in more detail in the bottom panels of 
Fig.m This break-up of the spectral weight into a quasiparticle peak along the needle (1) 
and a broad high intensity structure (2) is typical for the nodal ARPES measurements 
observed in a wide variety of cuprate compounds. As the break-up is produced by the 
charge-order gap, it should be considered as a direct manifestation of the local electron 
order in the cuprates. 

We note that the broad low-energy feature (marked as part (2) in Fig. [5]) extends 
in the fc-direction which is in contrast to the experimental observations. The reason for 
such a wide spread of this charge-order-induced feature (2) is a simplified approximation 
for the characteristic phonon frequency ujo = const. Within an Einstein approximation 
the electron-phonon contribution to the electronic self-energy S^f (a;) does not depend 
on the momentum vector k, which in turn leads to a wide spread of the self-energy part 
S2 and of the corresponding high spectral intensity region in the fc-space. To obtain 
better agreement with experiment, where the high-intensity structure is localized in k- 
space, one needs to consider a realistic description for the phonon dispersion which is 
beyond the scope of our current studies. 

Due to the doping dependence of ^2(^,0;), the high intensity features also depend 
strongly on doping which is demonstrated in Fig. [6l The high-intensity structure, at 
about —50 meV for x = 0.04, becomes smoother with increasing x up to x = 0.11 and 
in the overdoped regime {x = 0.2 in Fig. [6]). Moreover, as the local order is destroyed 
with increasing T, the contribution of ordered domains to I{k,uj) will decrease which 
can explain the fact of vanishing kinks for higher T discussed in Refs. |5] and |32] . 

The interpretation of ARPES intensities in terms of electronic inhomogeneities 
is a possible scenario not only for cuprates, but can be applied also for manganite 
compounds. In manganites, a significant electron-phonon coupling and an electronic 
phase separation of ferromagnetic metallic and charge ordered states is strongly 
supported by numerous experimental and theoretical studies [331 EH [35] . Consequently, 
the nodal ARPES features reported for Lai.2Sri.8Mn207 [36] can also be explained by 
the proposed approach for inhomogeneous states. 

In conclusion, we show that coupling to phonon modes leads to different spectral 
properties in ordered and disordered electronic states. It appears to be the key 
mechanism responsible for the main features detected in ARPES experiments of high-Tj, 
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cuprates. The broad character of the ARPES features can be related to the incoherent 
nature of the nodal non-Fermi-liquid state which forms essentially due to the electron- 
phonon coupling. In this context, the implications of electron-phonon coupling in charge 
ordered systems demonstrate in a novel way how collective modes can qualitatively 
change the fundamental properties of the electron liquid. 

This work was supported by Deutsche Forschungsgemeinschaft through SFB 484. 

References 

[1] Valla T et al. 1999 Science 285 2110 

[2] Lanzara A et al. 2001 Nature 412 510 

[3] Zhou XJ et al. 2005 Phys. Rev. Lett. 95 117001 

[4] Borisenko S Y et al. 2006 Pliys. Rev. Lett. 96 117004 

[5] Kordyuk A A et al. 2006 Phys. Rev. Lett. 97 117002 

[6] Damascem A, Hussain Z, and Shen Z H 2003 Rev. Mod. Phys. 75 473 

[7] Yunoki S, Dagotto E, and Sorella S 2005 Phys. Rev. Lett. 94 037001; Mishchenko A S and 

Nagaosa N 2004 Phys. Rev. Lett. 93 036402 
[8] Hengsberger M et al. 1999 Phys. Rev. Lett. 83 592 

[9] Valla T, Fedorov A V, Johnson P D, and Hulbert S L 1999 Phys. Rev. Lett. 83 2085 

[10] Byczuk K et al. 2006 Nature Physics 3 168 

[11] Pavlenko N and Kopp T 2006 Phys. Rev. Lett. 97 187001 

[12] Brown S E, Fradkin E, and Kivelson S A 2005 Phys. Rev. B 71 224512 

[13] Howald C et al. 2003 Phys. Rev. B 67 014533 

[14] Kivelson S A et al. 2003 Rev. Mod. Phys. 75 1201 

[15] Vershinin V et al. 2004 Science 303 1995 

[16] Kohsaka Y et al. 2004 Phys. Rev. Lett. 93 097004 

[17] Giustino L, Cohen M L, and Louie S G 2008 Nature 24 975 

[18] Held R, Bohnen K P, Zeyher R, and Manske D 2008 Phys. Rev. Lett. 100 137001 

[19] Cuk T et al. 2005 phys. stat. sol. (b) 242 11 

[20] Kohsaka Y et al. 2007 Science 315 1380 

[21] Tranquada J M et al. 1995 Nature 375 561 

[22] Hoffman 3 E et al. 2002 Science 295 466 

[23] Hanaguri T et al. 2004 Nature 430 1001 

[24] Lee J et al. 2006 Nature 442 546 

[25] Pan S R et al. 2001 Nature 413 282 

[26] Singer P M, Hunt AW, and Imai T 2002 Phys. Rev. Lett. 88 047602 

[27] Ofer R, Levy S, Kanigel A, and Keren A 2006 Phys. Rev. B 73 012503 

[28] Abbamonte P et al. 2005 Nature Physics 1 155 

[29] Eschrig M and Norman M R 2003 Phys. Rev. B 67 144503 

[30] Plate M et al. 2005 Phys. Rev. Lett. 95 077001 

[31] Abrikosov A A, Gorkov L P, and Dzyaloshinski I E 1975 Methods of quantum field theory in 

statistical physics (New York: Dover) 
[32] Johnson P D et al. 2001 Phys. Rev. Lett. 87 177007 
[33] Roder H, Zhang J, and Bishop A R 1996 Phys. Rev. Lett. 76 1356 
[34] Moreo A, Yunoki S, and Dagotto E 1999 Science 283 2034 
[35] Dagotto E, Burgy J, and Moreo A 2003 Solid State Communications 126 9 
[36] Mannella N et al. 2005 Nature 438 474 



